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Abstract 

Coexistence of various ordered chaotic states in a Hamiltonian system is studied 
with the use of a symplectic coupled map lattice. Besides the clustered states for 
the attractive interaction, a novel chaotic ordered state is found for a system with 
repulsive interaction, characterized by a dispersed state of particles. The dispersed 
and clustered states form an onion-like structure in phase space. The degree of 
order increases towards the center of the onion, while chaos is enhanced at the 
edge between ordered and random chaotic states. For a longer time scale, orbits 
itinerate over ordered and random states. The existence of these ordered states 
leads to anomalous long-time correlation for many quantifiers such as the global 
diffusion. 
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1 Introduction 

The ordering phenomena in Hamihonian systems is an interesting and important topic, 
for understanding the origin of chaos and order in nature. A traditional approach in 
the study of dynamical properties of Hamiltonian systems deals with invariant measure 
and ergodic properties. Here we focus on non-ergodic and non-stationary aspects of 
the systems and introduce various interesting phenomena. In a previous paper ([1]), we 
have shown that clustered motion is formed even starting from a random phase motion, 
in a class of Hamiltonian systems. In [1] it is shown that the fully developed chaotic 
( "random" ) state and clustered motion coexist in the phase space. The clustered motion 
is also chaotic with a different nature from fully developed ("random") chaos without 
structure. In the present paper we further explore various forms of order and chaos 
in a high- dimensional Hamiltonian system. In particular, we investigate more precisely 
properties of the clustered state and show that the clustered state is not uniform but 
actually consists of many clustered states with different degrees of order. We also show 
that, for repulsive systems, a new ordered chaotic state (called "dispersed order") exists. 
It will be shown that these clustered states, as well as the dispersed states, form an onion- 
like structure in the phase space. Following paper [1], we use a symplectic map system, 
since it is numerically efficient. Assume that we have particles on a unit circle and 
that the state of each particle is defined by its phase (position) 27r ■x(i) and its conjugate 
momentum p{i). Introducing the Hamiltonian 

Af TV- r ^ n=oo 

H = Y.P^?~ ,^ Ecos(27r(x(.)-x(j)) E Kt-n) (1) 

i=i (2vr)VA - 1 [^yi J „=-oo 

The temporal evolution of our model is given by; 

{xn{i),Pn{i)) ^ (^2;„+i(i),p„+i(i)),ii = 1,2,- ■■,A^, 
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K ^ 

Pn+l{i) = Pn{i) + - J=== sin 271 {x nU) - Xn{i)), , (2) 

zvrvJV — i 

where n represents the time step and i represents the particle index. When > 0, the 
interaction term 2itJn-i sin27r(x(j) — x{i)) between two particles i and j is attractive, 
while it is repulsive for i^' < 0. Total momentum J2f=i p{j) is a constant of motion of the 
model(0), so that the linear motion of the center of mass is separated from other dynamics. 
Thus we assume, throughout the paper, that the initial momenta satisfy J2jPn{j) = 
without any loss of generality. The temporal evolution rule is a canonical transforma- 
tion and the symplectic form is conserved: Thus 2N dimensional volume elements are 
conserved; 

d'^Pn A d^Xn = d^Pn+1 A d^Xn+i . 

As is shown in [|I| , the model shows clustering of particles for i^' > (attractive) case. 
Here we will also report a novel ordered state, for the repulsive case ( < 0). The order 
is characterized by a dispersed state of particles. Following [Q, we evaluate the degree of 
clustering and dispersion of particles with the following quantity; 



V def 
— 



N 




^exp (27rix„(j)) 



(3) 



■ ■ ■ if all x(j)'s are the same ( fully clustered) 

■ ■ - if x(j)'s are randomly distributed (4) 

■ ■ ■ if x(j)'s are evenly spaced 

The present paper is organized as follows; In §2, the repulsive case {K < 0) is studied. 
A novel order with dispersed state of particles is found with the use of the above order 
parameter Z. The onion structure of many ordered states is clarified in §3, with the 
use of Z. These ordered states have a finite lifetime, whose length increases as the orbit 
goes deeper into the onion structure. The ordered states are also temporally chaotic, 
whose strength varies by the degree of order, as will be studied in §4 with the use of 
Lyapunov spectra. The phase-space structure of many ordered states and the random 
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Figure 1: dispersed motion 



one are visualized in §5, by taking a small size system. For a much longer time scale, the 
orbit shows a switching process between ordered and random states. As will be shown in 
§6, this itinerancy over ordered states leads to anomalous behavior of many quantifiers; 
anomalous diffusion up to some time scale, stickiness at states without diffusion, as is 
shown in the diffusion coefficient and order parameter (Z) distributions. The residence 
time distribution of ordered states also exhibits a power-law behavior. Dependence of the 
order and randomness on the nonlinearity K is studied in §7 , where global diffusion of 
momentum is also discussed. Summary and discussions are given in §8. 

2 Repulsive case; Dispersed Order 

In the attractive case {K > 0), a clustered state is found when we start from an initial 
condition with a small momentum variation ( e.g., p{i) = random over [— pj„j/2,pj„j/2]; 
Pini is a small number , say O.2.), even if the initial phase x{i) is random (see [1]). In 
that case, most particles take nearby values, forming a cluster. In the present section we 
show there is a kind of ordered motion in the repulsive case {K < 0), corresponding to 
clustered state. Some timeseries for < are given in Fig.|l|, starting from an initial 
condition x{i) = random and p{i) = 0. Compared with the time series started from 
x{i) = random and p{i) = random, we see that these two time series have quite different 
character. A kind of order is seen in the former, but the form of the order is not so 
clear as in the attractive case just by the figure. Here Z provides a useful measure to 
see the existence of order. In the upper half of Fig.|l], the corresponding time series of Zn 
is plotted. The value is significantly less than unity (expected for a state with random 
phase). Through the repulsive interaction, particles push away each other. A possible 
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limiting case is 'equi-separation'; 



x{i) = const. + i/N 
p{i) = const. , 



(5) 



which is a stationary solution of our model and satisfies Z = . There are many 
other stationary states, e.g., x{i) = c for i:odd, x{i) = c + 1/2 for i : even, p{i) = 0, 



The above Z values, being much smaller than 1, suggest the existence of some order in 
the repulsive case, attained by the separation of particles. Here we call this state as a 
"dispersed" state . Thus the value Z gives an index of "order" in our model; 



Since clustered and dispersed states have many properties in common, we often call these 
states together as ordered states, here. On the other hand, a fully chaotic state with 
almost random phases { Z ^ 1) is called "random" state here, although its motion cannot 
be purely random, of course, due to the deterministic temporal evolution of our system. 
Simple linear stability analysis yields that the equi-separated stationary state (H) is stable 
for 



We have to note that ordered states are also chaotic, as will be discussed (see also [1]). 

3 Coexistence of Ordered states in Onion 

In , we treated the system as composed of two different chaotic seas; ordered and 
random. We see, however, the ordered state is actually composed of many different states 
depending on the initial conditions. To see this coexistence, we have used the following 




a state is 'clustered' if Z ^ 1. 



(6) 



while a state is called 'dispersed' if Z <^1. 
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Figure 2: Z vs pini 



Figure 3: Zn vs time 

one-parameter {pini) characterization of initial configuration; 

x{i) = random over [0, 1] and (8) 
p{i) = random over[-pj„j/2, pj„j/2]. 

By changing the initial momentum deviation pj„j, the nature of ordered states changes 
successively. In Fig. ^, the temporal average value of Z over initial 8000 steps is plotted 
with the change of initial momentum deviation p^j. We can see the coexistence of ordered 
states with different levels of order, coded by the value of Z. The ordered state exists 
only up to some value of p^j. This threshold and the degree of order, of course, depend 
on K as shown in Fig. ^, and will be discussed in §7. The average value of Z does not 
vary significantly within a family of (different random) initial conditions, if the parameter 
Pint is identical. Hence Pini is a relevant parameter to the degree of order in the system. 
To confirm that each ordered state is separated, we have also plotted the time series of Z 
in Fig. H, where Zn averaged over 4096 steps is successively plotted. At each time step, 
the value of Z is distinguishable by states. Thus many ordered states coexist in the phase 
space as distinct ones. In two examples in Fig. ^, the time series of Z switches from low 
values to those around 1.0, after a large number of time steps. This switch is a temporal 
transition from ordered states to the random state. As is discussed in [1], all parts in 
chaotic sea in the phase space are connected ^ , including all ordered states and the 
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Figure 4: Zn vs time (attractive) 



random state. Thus the ordered states are expected to have finite hfetime. Through time 
evolution an ordered state switches to the random chaotic state distinguished by the Z 
value around 1.0. ( If we wait for a long time, the reverse process is also possible, as 
will be discussed in §6.) In the attractive case, similar separation of the value Z is also 
seen in corresponding plots as Fig. ^. Howeever, the time course of switching process 
has a different character between attractive and repulsive cases. In the repulsive case, 
the switching process progresses rather rapidly, once it sets in. When the switch starts 
from an ordered state with a very low Z value, for example, the orbit passes through 
states with higher Z successively in a short time scale, till it reaches the random state 
(see Fig. |^). In the attractive case, on the other hand, the switching process is rather 
gradual, as is shown in Fig. ^. Crossover between ordered and random chaos is a novel 
type of chaos-chaos transition. As the initial randomness pmi decreases, the duration of 
such ordered state increases rather rapidly. The increase here might be estimated in a 
similar manner as given in the Nekhoroshev type argument [0, |19|, where the transition 
there is from an apparently regular state near a torus, to a fully chaotic state, while ours 
is chaos-chaos. The lifetime of an ordered state depends on its initial condition, which 
is well represented by the parameter mentioned above. We have plotted the lifetime 
with the change of Pmi in Fig. ^. In the repulsive case, the average transient time is 
roughly fitted with exp {const. / y/pmi) , down to some value of pini- (We have no reason 
to expect the divergence at Pini = 0, unless we take a regular initial condition for x{i), 
for example x{i) = i/N or x{i) = const.). This increase of lifetime is seen clearly in 
the repulsive case, while the increase with pini in the attractive case is saturated rather 
rapidly, and it is not easy to see a simple fitting form. The early saturation (without 
a simple fitting form) in the attractive case may be related with the gradual collapse of 
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Figure 5: Lifetime of clustered state 



ordered states. Summing up the section, ordered states coexist in the phase space, hke 
an "onion" structure. This type of onion structure may remind us of that supported by 
KAM tori and islands, as are often seen in pendula or the standard map 0. In contrast 
with these examples from low degrees of freedom, our "onion" consists of chaos, not of 
tori. 



4 Chaos in Ordered States 

As is noted in [1], the dynamics of the ordered state is also chaotic. Here we study how the 
strength of chaos varies with the degree of order in the clustered/dispersed states. For this 
purpose, we use the Lyapunov spectrum, a characteristic of asymptotic orbital instability 
of dynamical systems. It is a set of real numbers with 2N elements {Ai, ■ ■ ■ , A2Ar} and 
defined from an eigenvalue spectrum of the squared Jacobi matrix; 

J(.^ dof d{p{t),x(t)) 

{e^^i*, ■ ■ ■ , e^^^'^*} = eigenvalue spectrum of \j{t)J{t) as t ^ oo. (10) 

Note that, with this definition, a Lyapunov spectrum depends on the initial condition. For 
actual computation of exponents we use the standard method ||21|, ^ with Gram-Schmidt 



orthonormalization. We arrange the exponents as decreasing order Ai > A2 > ■ ■ ■ > \2N 
. Since A2Ar+i-i = — Aj due to the symplectic condition, only the bigger half of the whole 
spectrum is necessary. Some examples of Lyapunov spectra are plotted in Fig. ||, where 
initial conditions are taken as eq.(H). We note that all the exponents are much lower for 
Pini = than other cases. As Pini is increased, the exponents are shifted upwards, up to 
some value of and then decrease again to approach the spectra for the random state. 
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Figure 6: Lyapunov spectrum 



Figure 7: Max Lyap KS Z vs pini 



To see this pj„j-dependence clearly we have plotted the maximal Lyapunov exponent 
and Kolmogorov- Sinai (KS) entropy (estimated by Z^^i ^j) with the change of Pini in 
Fig. 1^ The enhancement of chaos at the border between random and ordered chaos is 
found both in clustered and dispersed cases. In the course of collapse of clustered states, 
chaos is enhanced as is shown in Fig. ^, where we have computed the local Lyapunov 
spectrum and Z, averaged over a given (finite) time steps (=4096 in the figure). From the 
computation, we have obtained the local Lyapunov exponents at each time step, as the 
short-time average of the exponents. As is shown in Fig. H, both the maximal Lyapunov 
exponent and KS entropy are gradually enhanced as Z is decreased ( that is, as the orbit 
goes to the outer part of the onion structure). With time increase, the local Lyapunov 
exponent and KS entropy take maximal values just before Z approaches unity. Reaching 
the random chaotic state, the Lyapunov exponents are again decreased to settle down to 
a constant value. The intermediate enhancement at the collapse is not clearly seen in 
the repulsive case, since the orbit passes through the outer part of onion within a short 
time scale. Chaos in a kicked system as ours often exhibits the diffusion in momentum 
space, as has been intensively studied in the standard map. Indeed the random state in 
our system shows diffusive motion, as is studied in |T^ and in §6. Is this true for chaos 
in ordered states? For the study of local diffusion in phase space, it is often useful to 
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Figure 8: Max Lyap vs time 



Figure 9: D(t) for clustered state 



introduce the following local diffusion coefficient defined by 



m ~ - Pn{^)f », (11) 



where the bracket << ■ ■ ■ >> represents the long-time average [II3, i-e. average over 
many ra's. If the diffusion in the phase space is normal, there exists a constant D^o = 
lim^^oo D(t) > 0. If the diffusion is fractional, D{t) = t~^, with some positive exponent 
6 < 1, which characterizes the stickiness of such diffusion. If the orbit is localized without 
any global diffusion, D{t) oc 1/t. For the ordered state, D{t) oc 1/t over all time steps 
within its lifetime (see Fig. This 1/t decrease is in contrast with the diffusion in the 
random chaotic state ( plotted in Fig. ^ for reference, by taking the initial condition 
Pini = !)• Thus there is no global diffusion in the momentum space for the ordered state, 
although the motion is chaotic. Indeed, the localization of orbits is necessary to have 
the onion structure, where the value of Z remains distinct by initial conditions. For the 
attractive case, D{t) for the clustered state is larger than that for the random state, if t is 
small. The momentum oscillates with a larger amplitude for the clustered state, but does 
not show global diffusion. On the other hand, D{t) ( even for small t) is much smaller in 
the dispersed state. It monotonically increases with p^i up to the random state. 
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5 Structure in the Phase Space 



To understand the onion structure of chaos and order, it is essential to explore the struc- 
ture in the phase space in detail. Since the "anatomy" of high-dimensional phase space 
is difficult to visualize, we study a rather low- dimensional case here, by restricting the 
number of particles to 4. (Due to the conservation law, the phase space is 3x2 dimension). 
From the 6 dimensional phase space, we take a two-dimensional slice for visualization. 
Here we sample the dynamics of initial conditions of 512 x 512 points in (p(l),p(2))-space 
by fixing x{l),x{2),x{3) and p(3). In addition, x(4) and p(4) are determined from the 
constraints J2'j=iP{j) = and J2'j=ix{.j) = 0. In Fig. Z values are plotted, obtained 
from the average over 256 steps starting from the corresponding phase space points. Cor- 
responding plots of the maximal Lyapunov exponent averaged over 256 time steps are 
given in Fig. \n\ . 

For the attractive case, the core of the onion is seen around {p{l),p{2)) ^ (0,0). 
The value Z increases monotonically as the point approaches the origin. We can see a 
threshold for the momenta below which the ordered states exist. Besides this expected 
structure, we have also seen regions corresponding to ordered states near a two-clustered 
state with 1:3 (around (p(l) = p(2) = p(4)) ^ (±1/3, ±1/3)). From the edge of the onion, 
broad resonant structures are emitted where Z values are rather low, and the Lyapunov 
exponent is rather high. This structure supports the enhancement of chaos at the edge 
between order and chaos. 

For the repulsive case, we have again observed the onion structures around {p{l) , p(2)) ^ 
(0,0), and the satellite structure (around (p(l),p(2)) ^ (±1/3, ±1/3)). The onion struc- 
ture, however, is not smoothly constructed as in the attractive case. The Z value changes 
in the momentum space not gradually, but it sensitively depends on the initial condi- 
tion even around (p(l),p(2)) ^ (0,0) ("chopped onion" structure). This may be the 
reason why the switching from dispersed to random states are rapid in the repulsive 
case. Furthermore some resonant structures are clearly visible satisfying the conditions 
p{l) = ±p(2), p(l) = ±2p(2), p(2) = ±2p(l), and so on. ( see also |T3| for the resonant 
structure in higher dimensional phase space, as well as the pioneering work ||31|| ). 
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Figure 10: Phase Space Slice with Z values 



Figure 11: Phase Space Slice with Lyapunov exponent 

A problem in the anatomy in the present section is that the size {N — 4) may be 
too small. This problem seems to be more serious in the repulsive case, since particles 
can often split into two groups with two elements or with one and three elements; such 
separation can lead to a a larger value of Z than 1, and bring about some difficulty in the 
distinction of ordered {Z < 1) and random states. The problem for the small size is also 
seen in the chopped onion structure, which is in apparent contradiction with the smooth 
onion picture in §3. We believe that the smoothness is attained with the increase of size ( 
even in the repulsive case) , where many degrees of freedom may smear out fine structures 
in the phase space. 

6 Order within Chaos 

The existence of ordered states affects the long time behavior of the (random) chaotic 
motion. The orbit visits many ordered states during the long term evolution. Switching 
between the random and ordered states can occur. How do the remaining ordered states 
affect the long-term statistical behavior of the dynamics? To address this question, we 
have studied the behavior of (1) local diffusion coefficient D(t), (2) local diffusion distri- 
bution, (3) residence time distribution at ordered state, and (4) distribution of local order 
parameter (Z). (1) Local diffusion coefficient D{t) D{t) is a good characteristic to 
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Figure 12: D(t) for random state 



search for order within chaos. As shown in Fig. |T2|, D{t) shows an anomalous power law 
decay D(t) oc up to some crossover time. This remnance of sticky behavior up to a 
large time suggests the long-time residence of the orbit near an ordered state. 

(2) Local diffusion distribution As another direct way to see the sign of ordered 
states, we have measured the distribution P{d) of short time diffusion coefficient for each 
particle; di(t) = {pn+tii) — Vnif))"^ ■ The distribution P{d) is given in Fig. |1^ for ordered 
and random states. The distribution is fitted in the following form; 

Pid) ^ / ^"^^^^ ^ (12) 

1 exp(— const. c?) for large d. 

For large \K\ (e.g., li^l > 0.7), the exponent a agrees with 1/2. The value 1/2 is easily 
explained by the central limit theorem: If we assume the position x{i) of particles are 
independent random numbers, we can expect that the distribution of the force term 
SjLi sin(27r(a;„(j) — x„(z)) obeys the Gaussian distribution. The (average of the) square of 
this variable gives the local diffusion di{t). If a variable z obeys the Gaussian distribution, 
the variable y = obeys the distribution ?/~^''^ exp(— const x y), thus leading to the 
above behavior. On the other hand, the exponent a is rather close to 1 for small K where 
ordered states exist in the phase space. The value a ~ 1(> 1/2) suggests the stickiness 
to states near D {). As is shown in the previous section, there is no diffusion in the 
phase space for the ordered states. Thus the enhancement of the exponent a means that 
the orbits are frequently stuck to ordered states. (3) Residence time distribution 
To see the residence at ordered states more explicitly, we have measured the residence 
time distribution at ordered states. Take a long time series and decompose the time 
interval into segments of 'clustered parts' and 'non-clustered parts' by properly chosing a 
threshold value Z , say, 1.05. ( We study the attractive case here). The length of each 
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Figure 13: P(d) Distribution of diffusion 



Figure 14: P(t) Residence time distribution 
segment represents a residence time at the clustered or non-clustered state. From the set 



of segments we can see the distribution of residence times as seen in Fig.^. Indeed the 
residence time distribution at ordered states obeys the power law distribution t~'^, c ~ 3.4, 
which means that no characteristic time scale exists in this region. For the random state 
we have exp(— c't) distribution, implying uncorrelated motion of particles. For the ordered 
state the exponent c is greater than 1 and the average residence time is finite. Thus the 
temporal switching between the ordered and random states continues forever. The above 
difference between the two distributions means that the switch between the two states 
is quite asymmetric. From random to ordered states the switch occurs through a kind 
of random trap to holes, while the reverse switch is rather gradual. The orbits slowly 
departs from the ordered states. This type of asymmetric switch is commonly observed 



in the chaotic itinerancy [27, 28 . (4) Distribution of local order parameter (Z) 
Another characteristic is a distribution of the order parameter Z over the phase space, 
where Z is calculated as finite time average. We have measured the distribution of Z 



averaged over a given time step r in Fig. 15. The distribution has a peak around 
Z ^ 1, with exponential tails to both sides. The tail extends deeper to Z < 1 (for the 
repulsive case) oi Z > 1 (for the attractive case), respectively. The inner part of "onion" 
structure is less frequently visited. From the distribution of Z, we have also studied the 
decrease of its variance with the increase of sampling time r. It is found that the variance 
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Figure 15: P(Z) Distribution of Z 



decreases with the coarse-grained time r as Var{Z) ^ r^^ for large K, in accordance 
with the central limit theorem, while it decays as r^^/^ for small K. This slow decay 
again suggests the sticky motion to ordered states. The scaling here corresponds to large 
deviation analysis The power law distribution in the residence time often leads 

to the anomalous scaling in the large deviation. Indeed Kikuchi and Aizawa [^] have 
shown a relationship between the power of residence time and the scaling for the variance, 
by taking a two-state semi-Markovian process. According to our numerical simulation, 
however, the relationship between c in (3) and the exponent here does not quantitatively 
agree with their theory. Following their theory, the central limit theorem with usual 
scaling ( oc r^^) is valid if c > 2 (as in our case), in contrast with our anomalous scaling 
with r. We believe that this discrepancy comes from the fact our ordered states are 
not single. An infinite-state ( rather than two-state) semi-Markovian process must be 
necessary for a quantitative analysis. Results of the present section are summarized in 



Table I. Table I 



Quantifiers 


Anomalous behavior 


cf. 


m 


crossover from to D^o 


Doo for h. n. 


local diffusion 






distrb. P{d) 




for h. n. 


residence time 




exp(— t) for 


distribution P(t) 


ordered states 


random states 


Var{T); 

variance of local Z 


r-V2 


for h. n. 



high nonlinearity ( large K) 



h.n. 
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Figure 16: Z vs K varying pini 



Figure 17: D vs K 

7 Parameter Dependence and global diffusion 

As our model (0) is integrable for K = 0, the parameter K can be regarded as the 
magnitude of perturbation to integrable model. As the perturbation K is changed, the 
degree of order in the ordered state varies, as well as the volume of the phase space 



supporting the ordered states. In Fig. |I6|, we have plotted the average value of Z over 
initial 8000 steps, starting from the initial conditions pmi = 0, .1, .2, ■ ■ ■ , 1.1. As is shown 
in Fig. |l^, ordered states disappear as \K\ gets larger. We note that the mechanism of 
disappearance is different between attractive and repulsive cases. In the attractive case, 
the value Z for clustered states decreases towards 1, the value for the random state, as 
K is increased (see Fig. 0). In other words, the "order" in a clustered state decreases 
with K till the state is absorbed into the random chaotic state. In the repulsive case. 



there remains a large gap in Z between the dispersed and random states (see Fig. jTB]). 
Thus the ordered states still remain as a structure even if \K\ is increased. Instead, the 
lifetime for these dispersed states decreases with \K\, till it is too short for the states to 
be observed as a temporally stable one, for large \K\. Of course the long-term behavior of 
chaos varies with the nonlinearity K. The diffusion constant D, estimated as limj^oo D{t) 
is plotted in Fig. O. First we have to recall that the coupling constant K is scaled by 



yiV— T so that the model is expected to show extensive behavior in a strongly chaotic 
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regime K^l. In the strongly chaotic regime, correlation among particles is negligible. 
Thus the force terms {2n^/N - 1)-^K J2f=i sin {27r(a;t(j) - Xt{i)))} can be approximated 
by stochastic variables independent of the system size A^. This approximation leads to the 
proportionality of diffusion coefficient to K"^, which is numerically confirmed for fi' ^ 1 
Due to the scale factor of y/N — 1 in front of K, the diffusion constant there approaches 
a size- independent value as is increased, (see Fig. In a smaller K regime, the 

diffusion constant D increases as K"^ with the exponent a different from 2. This power 



slowly decreases with the size A^, as shown in Fig. The estimated exponent a is given 
in Table II. We note here that this is not a Nekhoroshev-type dependence D oc exp{cK") 
1^, which is expected to hold for nearly integrable systems. The fractional power law 
dependence, which is also numerically observed in another model , reminds one of "Fast 
Arnold Diffusion" for intermediate perturbation strength ( see Chirikov and Vecheslavov 



25|). As the system size gets increased, however, the exponent a here is closer to 2, 
the value expected by random phase approximation, whereas the Fast Arnold Diffusion 
theory predicts a ~ 6.6. It may be necessary to study a regime with smaller K for such 



a large system (A^ > 40). Table II : D oc K"" 



Size A^ 


10 40 80 160 


Exponent a 


4.1 2.75 2.44 2.34 



8 Summary and Discussions 

In this paper we have discovered a new type of ordered states in Hamiltonian systems. 
This ordered state, termed as "dispersed order" , appears in a system with repulsive inter- 
action. It is sustained dynamically as in the clustered state previously discovered. In the 
dispersed order, particles are scattered in a well-organized manner. The order is charac- 
terized by the parameter = (1/A^)| Z]j exp(27rix„(j))p. For dispersed ordered states, 
the average of Zn takes much smaller values than unity. Discovery of repulsive order in 
our state reminds us of the Alder transition or Wigner lattice, where spatial order ap- 
pears through repulsive interaction between particles. In contrast with these established 
examples, our "dispersed order" state in the repulsive model is more subtle, neither with 
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a clear periodic structure in space, nor with a static order. The "dispersed order" does 
not form a regular lattice as in the Alder transition but rather resembles a liquid state. 
Although our example here has a global interaction among particles, it is rather straight- 
forward to introduce a model with a short-ranged interaction, which shows a repulsive 
order as in the present example. These examples suggest that the order formation in a 
Hamiltonian system with repulsive interaction is rather common in nature. The order 
in the repulsive interaction is also seen in dissipative systems. Indeed we have found a 
clustering state (with 3 clusters) sustained by the repulsion of each particle, for a glob- 
ally coupled (dissipative) circle map [^. Such a clustered state has a different nature 
from the attractive case. Both the clustered and dispersed states have many features in 
common as ordered states. These states in our Hamiltonian system form an onion-like 
structure in the phase space. The degree of order decreases as the initial momentum vari- 
ance is increased, till the "random" chaotic sea replaces the ordered state beyond some 
threshold for the momentum variance. The ordered states also show chaotic behavior, 
although the orbits do not show global diffusion in the phase space. This localization of 
orbits gives a basis for the onion structure. It is found that chaos is enhanced at the edge 
between the ordered and random states. Orbital instability is strongest at the transition 
from ordered to random states. Coexistence of ordered and random states in Hamiltonian 
systems is already reported in Here we have shown that the ordered state is actually 
a union of states with various degrees of order. The coexistence of many ordered states 
is important, particularly when we think of quantum versions of the models, since the 
wave function spreads over the various ordered and non-ordered states in the phase space. 
Although the clustered and dispersed states (for attractive and repulsive interactions re- 
spectively) have many chaotic features in common, there are some differences between the 
two, as is summarized in the following table. So far it is not clear how these differences 
are interrelated, and if they can be explained from the phase space structure. Table HI 
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Interaction 


Attractive [K > 0) 


Repulsive(ii' < 0) 


Z 


> 1 


< 1 


Collapse of order with time 


gradual 


sudden 


Behavior of life time 


early 


increase with 


with ^ 


suppression 


exp(-Pi™' ) 


Collapse with 


merge into 


lifetime 


the increase of -f^ 


random states 


goes to 


-D(t) for small t 


larger than 


smaller than 




the random state 


the random state 


Phase space structure 






for iV = 4 


Smooth Onion 


Non-smooth 



Since all chaotic 



orbits are connected for a Hamiltonian system with many degrees of freedom, an ordered 
state has a finite life time before it is switched into the random state. Conversely, an orbit 
in the random chaotic state visits ordered states. The existence of the ordered states bring 
about anomalous long-time behavior in dynamics. Local diffusion (up to some crossover 
time) and the residence time distribution at ordered states show anomalous power-law be- 
havior, as well as the distributions of local diffusion coefficients and of the order parameter 
Z. The above behaviors are the manifestation of the existence of long-time correlation. 
Such power-law correlation implies the 1//" spectra for the Fourier transformation of 
autocorrelation of dynamic variables {a = 2 — 6; see Jl^)- We note that these long-time 
correlations always appear when the coupling K is small. In other words, we can expect 
that long-time correlation generally appears in a Hamiltonian system ( which is also true 
of a continuous-time case; see the later argument). Thus we can expect that 1//" behav- 
ior generally appears in a high-dimensional Hamiltonian dynamical system. This explains 
at least some of the origins of 1//" spectrum in nature, in particular in the fluctuation 
around equilibrium states p6| . The exponent a indeed approaches 1 as the coupling K is 



decreased [|T^]. The lifetime of the clustered and dispersed states increases with li^l 0. 
Since a proper limit with \K\ — > gives a flow system of a time-independent Hamilto- 
nian, we can expect that clustered and dispersed order are more frequently observed in a 
Hamiltonian system with a continuous time. In the continuum limit, dependence in 
the present paper can be related to the energy dependence. In the limit, it is expected 
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that many ordered states ( in the onion structure) exist as distinct stable states, depend- 
ing on the energy of the system. In this case our resuhs imply that the degree of order 
decreases up to some energy, beyond which the order collapses. Switching among ordered 
states through high-dimensional chaotic states have been extensively studied in dissipa- 
tive systems, as chaotic itinerancy ||2^. Switching between our clustered motions through 
random chaos provides an example of chaotic itinerancy in Hamiltonian systems. Similar 
ordered motion is seen in molecular dynamics simulations for glass |^ and water ||30| . 
Both the dispersed and clustered states will hopefully be found in other physical systems, 
such as gravitational systems, microclusters of atoms, colloids, and so on. Phenomena 
which have been explained by stationary solutions so far may actually be ordered states 
sustained by chaotic motion as in our example. This paper is dedicated to the memory 



of Jeff Tennyson, one of the pioneers in chaos and diffusion in Hamiltonian systems with 
many degrees of freedom Unfortunately I (=KK) only had three chances for dis- 



cussions with him at Berkeley and Los Alamos; or I should thank that there were such 
chances at all. I was always impressed by his deep thought. Besides this scientific im- 
pression, I somehow felt that Jeff might have had some difficulties in adapting to his own 
society, and felt that he might have much in common with the Eastern way of thinking 
and living. We would like to thank Dr. Y.Aizawa, K. Ikeda, K. Shinjo, T.Yanagita, T. 
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Figure Captions 

Fig. 2 Typical examples of overlaid timeseries of our model (0), with corresponding time- 
series of Z„. System size = 16, i^' = —0.4. (a) dispersed motion : initial con- 
dition x{i) = random, = 0. Plotted over the time steps 600 - 1200. (b) 
dispersed motion : initial condition x{i) = random, = random over [-0.2,0.2]. 
Plotted over the time steps 600 - 1200. (c) random motion : initial condition 
x{i) = random, = random over [-0.8,0.8]. Plotted over the time steps 800 - 
1000. 

Fig. ^ Average value of Z plotted as a function of pmi- K = —.5, —.3, —.1, .1, .2 and .5. 

= 16. The value Z is obtained from the average of initial 8000 steps, starting 
from a randomly chosen initial condition parametrized by pini- 

Fig. ^ Time series of Zn, for different values of pmi- N = 16, and K = —A. Plotted per 
4096 time steps. Pini =0,.l.,2.,3.,and .4 (from bottom to top). 

Fig. ^ Time series of Z„, starting with a random initial condition of Pini = 0.1. K = .15, 
and A^ = 16. Plotted per 4096 time steps. 

Fig. ^ Semilog plot of the average lifetime of ordered state versus l/i/Pmi- The average 
is taken from 80 samples (for .5 < Pini < .8) with different random configuration 
with given pini, and from 40 samples for .2 < pini < .5 , and from 20 samples for 
.15 < Pini < .2. The lifetime is estimated as the time step at which Zn exceeds 0.99 
for the first time. K = —0.4 and A^ = 16. 

Fig. ^ Lyapunov spectra Aj for the model N = 16. Spectra are obtained from the aver- 
age over initial 10000 steps, (a) K = —.2. Spectra from pini = 0, .1, .3, .5, .6, .75, 1, 1.5 
and 2 are overlaid, (b) K = .2. Spectra from pini = 0, .1, .2, .3, .5, .6, .75, 1, 1.5 and 
2 are overlaid. 

Fig. ^ The maximal Lyapunov exponent ( thick line: multiplied by 10 times for the conve- 
nience of scale), KS entropy (dotted line) , and the average value of Z (broken line) 
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are plotted as a function of initial pini- N = IQ. Obtained from the average over 
initial 10000 steps. K = .2 { KS entropy is multiplied by 5 for scale). 

Fig. ^ Temporal evolution of local Lyapunov exponents. The maximal Lyapunov exponent 
( thick line: multiplied by 10 times for the convenience of scale), KS entropy (dotted 
line; multiplied by 5 times) , and the average value of Z (broken line) are plotted 
as a function of time per 4096 steps. All the values are the average for 4096 time 
steps at each time. K = 0.15, Pmi = 0.1, and = 16. 

Fig. ^ Local diffusion coefficient D{t) for ordered states; Obtained with 10 - 100 sequential 
sampling, starting from the initial condition with pini = 0, .1, .2 and 1. 



Fig. A 2 dimensional surface of section (p(l),p(2)) of the model (0). Number of particles 
is A^ = 4, and couphng constants are (a) K = —0.05, (b) K = 0.05, (c) K = 0.2. In 
the (4x2=) 8-dimensional phase space, the section is taken by setting the following 
6 constraints; (b) x(l) = -0.075 = -x(4), x{2) = -0.025 = -x(3), p{3) = 0, 
and p(4) = — p(l) — p(2) so that the center of mass is fixed; J^i^i^) = J2iP{i) = 
. The last two constraints come from conservation of the total momentum so that 
the points shown in the figure has the same value of total momentum, which is a 
conserved quantity. We set 256x256 (for K = 0.05) or 512x512 (for K = -0.05, 0.2) 
lattice points on the 2-dimensional section, which are taken as initial conditions for 
the time evolution. The value of Z averaged over 10000 steps are plotted with a 
gray scale corresponding to the initial condition of (p(l),p(2)). 



Fig. 11 Maximal Lyapunov exponents are plotted depending on the initial condition. The 



same slice as Fig. |T0| is used. Starting from the initial condition corresponding to 
the point in the slice, we have measured the maximal Lyapunov exponent averaged 
only over 256 time steps, and they are plotted with the gray scale. K = 0.2. 

Fig. |l^ Local diffusion coefficient for the random state; Obtained with the average over 200 
sequential sets of data ( totally from 200 x t steps), starting from a random initial 
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condition with pini = 1. N = 40. Log-Log plot. K = .4, .3, .2, .15 from top to 
bottom. 



Fig. Distribution of diffusion coefficient d{i) = — The average time r is 

chosen to be 512. = 16, and Pini = 1. Sampled over 10000 times per r steps. 
Log-log plots, (a) K = 0.15 (b) K = 0.8. 

Fig. |l^ Residence time distribution at ordered and random states (a) for ordered [Z > 1.05), 
(b) for random {Z < 1.05) states. = 8, and K = 0.2, Total length of time series 
= 10^. Note that the distributions (a) and (b) are independent each other. 



Fig. Distribution of Z averaged over r time steps. = 16, and pmi = 1. Sampled over 
10000 times per r steps (a) K = -0.4, with r = 2048 (b) K = 0.15 with r = 1024. 



Fig. The average value of Z over 8000 steps, with the change of K, starting from a 
random initial condition with pmi = 0, .1, .2, ■ ■ ■ , 1.1. = 16 



Fig. The diffusion constant D vs K. D is estimated by the asymptotic value of D{t) 
converging with the increase of t, from the plot of D{t) corresponding to Fig. |1^. 
Obtained with 200 sequential average, starting from a random initial condition with 
Pini = 1. N = 10,40,80 and 160. Log-Log plot. 
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